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^ ! Abstract 



OO 



We perform simultaneous calculations of the radiative decays of scalar 
^ ! mesons /o(980) 77, ao(980) 77, vector meson 0(1020) 7/o(980), 

7ao(980), 77r'', 77/, 77/' and tensor mesons 02(1320) — > 77, /2(1270) — > 77, 
/2(1525) 77 assuming all these states to be dominantly the qq ones. A 
00 , good description of the considered radiative decays is reached by using almost 

' the same radial wave functions for scalar and tensor mesons that supports the 

idea for the /o(980), ao(980) and 02(1320), /2(1270), /2(1525) to belong to 
^ , the same P-wave qq multiplet. 

ix. 

^ : 1 Introduction 
> 

' Despite a long history of the P-wave qq multiplet the problem of definition of qq 

;h \ scalars is still a subject of lively discussion, see e.g. [@, U, H and references therein. 

Radiative decays of mesons may serve as a useful tool for the study of qq structure of 
mesons, in particular, P-wave qq component in /o(980) and ao(980). In this way, it 
is rather important to investigate simultaneously the other mesons which belong to 
the P-wave qq multiplet, namely, tensor mesons a2(1320), /2(1270) and /2(1525). In 
the present paper, combined calculations of the decays ao(980) — >• 77, /o(980) — > 77, 
02(1320) — > 77, /2(1270) — 77, and /2(1525) 77 are carried out assuming the qq 
radial wave functions in these mesons to be nearly the same. 

Radiative decays of the 0-meson are another source of important information 
on scalar mesons. We have calculated the decay processes with the production of 
mesons belonging to scalar and pseudoscalar sectors: 0(1020) 7/o(980), 7ao(980) 
and 0(1020) — >• ■j7io,'yrj,'yri'. These latter, of the type of V ^ 7P, are the classical 
reactions, which had been used rather long ago for the determination of qq structure 
of vector and pseudoscalar mesons |^ . 
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We believe that simultaneous description of the processes S — >• 77, T — > 77, 

V — 76* and V — *■ 7P is a necessary test for the whole calculation procedure and 
determination of the P-wave qq multiplet. 

In calculations of the decay form factors we use spectral integration over qq states 
together with the light-cone wave functions for the qq mesons; the method of the 
spectral integration for the form factor amplitudes has been developed in a set of 
papers [§,0,1,0. 

In Section 2 we present necessary elements of the technique for the calculation 
of radiative decay amplitudes. The detailed presentation of the technique for the 
description of composite qq systems can be found in Refs. [0, ||, §], where the 
pion form factor was studied together with transition form factors 7r° 7(Q^)7, 

V ~^ 7(Q^)7 cind rj' 7(Q^)7- The method of spectral integration works for 
the form factor amplitudes which obey the requirement of analyticity, causality and 
gauge invariance. The used technique allows one to introduce the composite-particle 
wave functions and perform calculations in terms of the light-cone variables. 

It is worth noting that this calculation technique for the processes involving 
bound states has a broader applicability than for the qq systems only: in |^ this 
method got its approbation by describing the deuteron as a composite np system. 



then this very technique was applied to heavy mesons [0. For the convenience of 
a reader, in Section 2 we recall briefly the basic points of this approach. Then we 
give necessary formulae for the calculation of partial widths for the decays V 76", 
V — >• 7P, S* — >• 77 and T — >• 77. 

Results of the calculation are presented in Section 3. First, we discuss the de- 
cay 0(1020) — s> 7/0(980). Our calculations show that the data on branching ratio 
PP(0(98O) 7/o(980)) = (3.4 ± OAtH) ■ 10"^ |n|, [g may be described assuming 



the qq structure of /o(980) and varying ss and nn = (uu + dd)/\^ components 
in a broad interval. For the flavour wave function written as ip fiavour[fo{^^0)] = 
nncos(f + sssiny) the agreement with data is reached with 25° < \(f\ < 90°. 

In Section 3, we calculate also partial widths for the decays 0(1020) •yr], ■jrj', 77r°, 
with the same technique as has been used for the reaction 0(1020) 7/o(980) 
and with the same parametrization of the 0-meson wave function. The calcula- 
tions demonstrate a good agreement with data as well. It should be stressed that 
in fact the decays 0(1020) — > 7?7,7?7' are calculated without any free parameter: 
these decays are governed by ss components in rj and rj' which are well known; the 
wave functions of the basic pseudoscalar and vector mesons are also known, see e.g. 



0, 0. So, the calculations of the decays 0(1020) — *• 7//, 7//' are needed for the 
verification of the method only, and the results provide us with a strong argument 
that the applied method for calculation of radiative decays of qq mesons is wholly 
reliable. The decay 0(1020) 'yir^ allows us to estimate the admixture of the 
nn component in 0(1020). With the flavour wave function of the 0(1020) written 
as V'/iaDour (0(1020)) = SS COS (pv + nnsinipv, we have \(pv\ < 4°. Partial width of 
0(1020) 7ao(980) is also proportional to the probability of nn component in the 
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0(1020); we discuss this decay in Section 3 as well. 

Two-photon radiative decays provide important information about P-wave qq 
mesons; the technique for the calculation of scalar- and tensor-meson decays, — 77 
and T 77, is presented in Section 2. 

Under the assumption of the qq structure of ao(980) and /o(980), analysis of 
partial widths ao(980) — > 77 and /o(980) 77 has been performed in 0. The 
data for ao(980) — > 77 are in a reasonable agreement with calculation. Concerning 
the /o(980), the extraction of the signal /o(980) 77 from the measured spec- 
tra 77 nir faces strong interference "resonance + background", thus resulting 
in uncontrollable errors (see, for example, the K-matrix calculation of the S'-wave 



spectra 77 vrvr [0). The recently obtained partial width r(/o(980) 77) 
0.28loj3 keV [|T3| is by a factor 2 smaller than the averaged value reported previously 



(0.56 ± 0.11 keV [0). In Section 3 we re-analyse the decay /o(980) — » 77 using 



new data for the width. Thus we get two allowed intervals for the nn/ss mixing 
angle: 80° < ip < 93° and (—54°) < (f < (—42°). The restrictions for mixing angle 
ip which come from the combined analysis of radiative decays 0(1020) —>■ 7/o(980) 
and /o(980) — > 77 are discussed in Section 4, we have two solutions for ip: 

V? = -48°±6°, if = 86° ±3°. (1) 



For positive mixing angle, the analysis gives us strong restriction for the value of 



radius of the /o(980): -Rlcggo) ^ GeV ^ (remind, the pion radius squared is 



Rl ~ 10 GeV-2). 

In Section 3 we discuss also the results of the calculation of tensor-meson two- 
photon decays: 02(1320) 77, /2(1270) — ^ 77 and /2(1525) 77. The form 
factors of corresponding transitions depend strongly on the choice of the vertex 
T — i> qq. In line with the qq classification, we perform comparison with data for 
the vertex which is related to dominant qq P-wave. The results are in reasonable 
agreement with the measured partial width r(a2(1320) — 77), calculations being 
carried out with the wave function of 02 (1320) whose characteristics are close to 
those of ao(980), namely, ^02(1320) — -^ao(980) — 7 — 12 GeV~^. Description of the 
two-photon decays of /2(1270) and /2(1525) fixes nn/ss ratio for these states. With 
flavour wave functions written as V'/«ai)owr [/2(1270)] = nncosipT + sssimpT and 
Ip fiavourif 2(^525)] = —nnsimpT + sscosipT, simultaneous description of the data 
can be reached at Pj2(i270) — -^/2(i525) — 7 — 10 GeV~^ requiring either ipx — 0° or 
cpT ^ 25°. 

Simultaneous description of radiative decays of scalar and tensor mesons, /o(980), 
ao(980) and 02(1320), /2(1270), /2(1525), with the use of similar radial wave func- 
tions, argues in favour of their belonging to the same P-wave qq multiplet. 
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2 Radiative decays in the framework of spectral 
integration technique 



Given here are the formulae for partial widths of the radiative decays: V — > ^yS, 
V — *• ■jP, 5 — 77 and T —> 77. Using as an example the reaction V — > 75, 
we present necessary elements of the spectral integration technique applied for the 
calculation of transition form factors. 



2.1 Moment-operators for the transition amplitudes S 

77, T 77, V and V 7P 

For amplitudes under consideration, we present the moment operators for two- 
photon decays of the scalar and tensor mesons: /o(980), ao(980) 77 and /2(1270), 
02(1320), /2(1525) 77, and for radiative decays of the 0-meson: 0(1020) 77r°, 
7?7, 777', 700 (980), 7/0 (980). Systematic presentation of the moment operators is 
given in [|17[] . 



2.1.1 Transition amplitude 5 ^ 77 

The transition amplitude S — > 7^(5'^) 7^(5''^) for the transversely polarized photons 
reads: 

A^p = e'Fs^,,{q',q")gi^ . (2) 

Here e is the electron charge {e'^/Air = a = 1/137); the indices a, (3 refer to the 
photons; q and q' are the photon momenta. The metric tensor g^^ works in the 
space orthogonal to p = q + q' and q: 

±± _ q±aq±f3 PaPl3 _ ± ± _ PaPa' , . 

9al3 — 9olI3 2 2 ' ^-L" ~ 9aa'1a' , Qaa' ~ 9aa' o ' ^ 



2.1.2 Tensor meson decay amplitude T ^ 77 



The T ^ 77 decay amplitude has the following structure: 



A 



s'^'JUp, q) F^ZUO, 0) + s':'Jp, q) F^Z(0, 0) 



(2) 



.(2) 



(4) 



where 5", 



(0) 



}^u,ai3iPyl) S'^d S^J^^i^{p,q) are the moment operators, indices a,P refer to 
photons and /i, u to the tensor meson. Two transition form factors for the trans 
versely polarized photons T — >• 7_L(g^)7±(g'^), namely, F^°^^^(g^, g'^) and F^}^^^ 
depend on the photon momenta squared q^ and g'^; the limit q^ = g'^ = corre- 
sponds to the two-photon decay. 



{q\q'^) 
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The moment-operators read: 
and 

s'^^J, apip^ q) = attduft + at^hua - atudip ■ (6) 

The moment-operators are orthogonal in the space of photon polarizations: 

c(o) c(2) _ n 

2.1.3 Transition amplitude V 'jS 

Transition amplitude V l±{(l'^)S for the transversely polarized photon takes the 
form: 

A,^ = eFv^,s{q')9i^a ■ (7) 

The index a refers to the photon and fi to the vector meson; p and q are the momenta 
of the initial vector meson and photon. The metric tensor g^^ works in the space 
orthogonal to p and q. The limit q^ = corresponds to the radiative decay of vector 
meson. 



2.1.4 Transition amplitude V — > 7P 

The spin operator for the amplitude V — >• 7P contains antisymmetric tensor e^uap, 
and the amplitude has the following structure: 

A^lOi = e ^^iauiU2PviQu2Fv-*"/p{Q^) ■ (8) 
The notations are the same as for V 76*. 

2.2 Form factor for the radiative decay V ^ ^±{q^)S . 

The method of calculation of the three-point form factor amplitudes in terms of the 
spectral representations over the qq intermediate state masses was developed in 
Here we give schematic presentation of the method using as an example the reaction 

V ^ 7±(g')^. 



2.2.1 Double spectral representation of the form factor. 

Assuming the qq structure for the initial (V^) and final (S) mesons, the amplitude of 
the decay V is determined by the subprocesses V ^ qq and qq S, with the 

emission of 7(5'^), see Fig. la. Corresponding three-point loop diagram is calculated 
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using double spectral representation over intermediate qq states: they are marked 
by dashed lines in Fig. 16. 

To be illustrative, let us start with the three-point Feynman diagram. For the 
process of Fig. la one has: 

/I (Feynman) _ / ^ ^V~*"iS fiu ^ /q\ 

~ J z(27r)4 ^ (m2 - A;2)(m2 - k[^){m^ - k^) ^ " 



Here ki, k[, ^2 are quark momenta, m is the quark mass, and Gy, Gs are quark- 

7(99) 



meson vertices; the quark charges are included into Zy'^l g. The spin-dependent 



block reads: 

SZ^''^ = -SP liK + mhtih + mh^{-k + m)l , (10) 



where the Dirac matrices 7^ and 'j^ are orthogonal to the emitted momenta: J^q^ = 
and 'j^Pu = 0. 

To transform the Feynman integral (^) into double spectral integral over invariant 
qq masses squared, one should make the following steps: 

(i) consider the corresponding energy-off- shell diagram. Fig. 16, with P2 = (^^ + 
^2)^ > 4m2, = {k[ + k2y > -iini? and fixed momentum transfer squared q^ = 

(p - p')^ 

(ii) extract the invariant amplitude by separating spin operators, 

(iii) calculate the discontinuities of the invariant amplitude over intermediate qq 
states marked in Fig. 16 by dashed lines. 

The double discontinuity is the integrand of the spectral integral over P^ and 
P'2. Furthermore, we put the following notations: 

For the calculation of discontinuity, by cutting the Feynman diagram, the pole terms 
of the propagators are replaced with their residues: (m^ — k"^)^^ 5{Tm? — k"^). So, 
the particles in the intermediate states marked by dashed lines I and II in Fig. 16 are 
mass-on-shell, kf = k^ = kf = w?. As a result, the Feynman diagram integration 
turns into the integration over phase spaces of the cut states. Corresponding phase 
space for the three-point diagram reads: 

d^P, P'; fci, k2, k[) = d^P; A;i, k2)d^P'; k[, A;^)(27r)32A;2o5^'Hk2 - k2) , (12) 

where the invariant two-particle phase space (i$(P; ki, ^2) is determined as follows: 

The last step is to single out the invariant component from the spin factor (|10]). 
According to (|^, the spin factor ( pUD is proportional to the metric tensor, Sj^'^"'^^ ~ 
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df^vi which works in the space of the intermediate state momenta. Then the spin 



factor Sy^^g, determined as 



S%--'^=9t^Sy^,s{s,s',q'), (14) 

is equal to: 



4ss'q;(s, s', q^^ 
s + s' — 



S'v^^5(s, s', g ) = -2m |^4m +s-s+q , J — b~ ) ' ^^^"^ 

I I q^{s + s'-q^) 

a{s,s,q ) - 



2q^{s + s') - (s - s')2 - g4 

Recall that, when going from (p!OD to (|15]), we use the mass-on- shell relations {kik2) = 
s/2 - m^, (fc'i^a) = s72 - m^, and (fc^/ci) = - g72. 

The spectral integration is carried out over the energy squared of quarks in the 
intermediate states, s = P"^ = {ki + k2Y and s' = P'^ = {k[ + A;2)^, at fixed 
g2 — (^p' _ py^ spectral representation for the amplitude Ay-^-ysiq"^) reads: 

oo , oo 



Av^^q^) =1-1 -^iP^ X (16) 



(h 7 ds' Gv{s) Gs{s') 
n J 

X j d^P, P'; k,, k2), k[ 4%s(s, s\ q')ZfXs ■ 

The spectral representation of the amplitude Ay^^siff) gives us the invariant part 
of A^^^^'^™'^^\ Eq. (H), when the vertices, Gy(s) and Gs{s')^ are constant. Generally, 
the energy- dependent vertices can be incorporated into spectral integrals. According 
to [^], the form factor of a composite system can be obtained by considering the 
two-particle partial- wave scattering amplitude 1-1-2 — 1 -|- 2: the pole sigularity 
of this amplitude corresponds to the composite system. The amplitude for the 
emission of a photon by the two-particle-interaction system has two poles related to 
the states "before" and "after" electromagnetic interaction, and the two-pole residue 
of this amplitude provides us the form factor of composite system. When a partial- 
wave scattering amplitude is treated using the dispersion relation N/ D-method, the 
vertex G{s) is determined by the A^-function: the vertex as well as A^-function have 
left-hand-side singularities which are determined by forces between the particles 1 
and 2. 

It is reasonable to name the ratios Gv{s)/{s — rn?) and Gs{s')/{s' — rn?) the 
wave functions of vector and scalar particle, respectively: 

Gv{s) Gs{s') 

2=^v{.s), - = ^s{.s). (17) 

s — s — m"^ 

Working with Eq. (|T6|) one can express it in terms of the light-cone variables. 
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2.2.2 Light-cone variables 

One can transform Eq. ([I6|) to the hght-cone variables, using the boost along the 
2;-axis. Let us use the frame where initial vector meson moves along the z-axis with 
the momentum p ^ oo: 

P=(p+f,0,p), P' = (p + l±^,-q^,p). (18) 
2p 2p 

In this frame the two-particle phase space is equal to 



IQn^ X1X2 



X Js-^^^^±^-^^^^±^] , (19) 

\ Xi X2 J 

where Xi = ki^/p, and the phase space for the triangle diagram reads: 

d^P, P'; ki, k2, k[) = ^^^^^d^k,^d'k2±5il - xi - X2)5(2)(kix + k2^) 

167r X1X2 

^ / _ + kl^ _ m'^ + klA ^ / , ^ 2 _ + (k^ - q^)^ _ m^ + fc|^ \ 
V Xi X2 J \ Xi X2 

Then the amplitude V —>■ 'y{q'^)S is written as 



(20) 



vwy; r dr r 
Av^,^,^)siq') = -tFt / -n ^ / d'k^Ms)Ms')Sv^ysis, s\ q') , (21) 



7 ill) 

IGtt'^ J xil — x)2 ^ 
^ ' 

where a; = ^22/^ , k^ = k2_L, and the qq invariant masses squared are 

_ + k'i , _ + (kx + xcii_f 

x(l — a;) x(l — x) 

2.2.3 Charge factors 

Charge factors for the nn and ss components in the transition — >■ 7/0 are deter- 
mined as: 

r^{nn) _ r)Ann) _ 1 



7(ss) _ r,>(ss) _ o„ _ _ 
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2 
3 



where Cjp^lfg is the following convolution: C^'^T-^/o = (uu+dd) / V^-eq-{uu+dd) / \p2 = 
i^u + ed)/2. Here e„ and are charges of and d quarks. The factor 2 in ( pSj) is 
related to two possibilities for photon emission, namely, from quark and antiquark. 
Likewise, for the process — > 709, one has: 
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2.2.4 Meson wave functions 



To calculate the form factors, one should define meson wave functions. The simplest 
parameterization is exponential one: 



(25) 



The parameters by and 65 characterize the size of the system, they are related to the 
mean radii squared, Ry and of the mesons. At fixed Ry and i?| the constants 
Cy and Cs are determined by the wave function normalization, which itself is given 
by meson form factor in the external field, -Fmeson(Q'^), and at small the form factor 
is: 
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:R„ 



(26) 



The requirement -Fmeson(O) = 1 fixes the constant Cmeson in (P^D, while the value 

^meson IS dlrCCtly related to &meson- 

In terms of the light-cone variables, the form factor -Fmoson(Q'^) reads: 



dx 



IQtt^ J x(l — xY 

^ ^ 

where 5'^eson is determined by the following traces: 

2P^SPis, s', q') = -Sp [ik[ + mh^Ch + ^)(-^2 + m) 
2P^S^''\s, s', q') = -isp [f',^{k[ + mhiih + mh^{-k + m) 

and the orthogonal components entering (28) are as follows: 



(Pg) 



7a 



7a Pa p2 
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P' 

'y - P' — 

la ^ a p,2 



(27) 



(2^ 



(29) 



where q = k[ — ki. When determining Sy^\s , s' , q'^) , we have averaged over three 
polarizations of vector meson that results in the factor 1/3. 

The functions Sg^\s, s', q^) and Sy^\s, s', q^) are equal to: 



Sf\s, s', q') = a{s, s', q^){s + s' - 8m' - q') + q' , 
St^'\s, s', q') = \ [a(., s\ q')is + s' + Am' - q') + q\ 



(30) 



where a(s, s', q ) is given in (p!5|). 
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2.2.5 Partial width 



The decay partial width is determined as 

mvTv^ys = \J d'^iPV, q,p's)\^f^u\'^ ■ (31) 

Here the averaging over spin projections of the meson and summing over photon 
ones is carried out (summation over photon spin variables results in the metric tensor 
g^^); the two-particle phase space for the radiative decay V ^ ■y + S is equal to: 
/ d^{pv', q,p's) = {fny ~ 'n^|)/(167rmy-). Partial width in terms of the form factor 
reads: 

mvTv^^s = -« '"^~"'^ |Fy^^5(0)P . (32) 



2.3 The process V ^ 

The light-cone representation of the transition form factor V —>■ 'y±{q'^)P reads: 



loTr"* J x 1 — x r J 
^ ^ 

The spin factor for pseudoscalar mesons is determined by 

-Sp i'y5{k[ + m)'y^-^{ki+m)'y'^{-k2 + m) = e^au^uiPu^quiSv^-fpis, s' , q^) (34) 

where q = P — P' and qala'^ = 0, P'^l'^ = 0. The spin factor is equal to 

S'y^^p(s, s', q^) = Am . (35) 

Charge factors for the considered radiative decays are as follows: for the ss com- 
ponent in the reactions — > 'yri,'yri', Z^^^^ = Z^^^^, = —2/3, and for the 7r° and 
ao(980) productions, Z^^^^o = Z^^.,ao{980) = 1- 



Partial width for the decay V ■yP is equal to 

1 /■ ,^ . , X , . ,o 1 ?7T.?/ — rrii 



mvTv^-yP = - / d<^{pv;q,Pp)\A^u\^ = -a ^ , ^ |Fy^^p(0)p . (36) 



2.4 Processes 5 ^ 77 

Our calculation of the two-photon decays of scalar and tensor mesons is based on 
the method developed in |^ for the study of the pseudoscalar meson transitions 
7r° — > 7(5^)7, r] 7(^7^)7 and rj' l{q^)li and for the photon we use quark- 
antiquark wave function which was found in We perform the calculation of the 
scalar and tensor meson transition form factors V l{q^)l and T — 7(^^)7 in the 
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region of small g^; these form factors, in the limit — 0, determine partial widths 
5* — >• 77 and T — > 77. 

The transition form factor qq — meson —>■ 7(^^)7 is determined by the three-point 
quark loop diagram of Fig. 16 type that is a convolution of the gg-meson and photon 
wave functions, ipqg ® ip-y. Following we represent the photon wave function as 
a sum of two components which describe the prompt production of qq pair at large 
s' (with a point-like vertex for the transition 7 — > gg) as well as at low-s' region 
where the vertex 7 — > gg has a nontrivial structure due to soft gg interactions. The 
process of Fig. 16 at moderately small |g^| is mainly saturated by the contribution 
of the low-s' region, in other words, by soft component of the photon wave function. 
The soft component of the photon wave function was restored in on the basis of 
experimental data for the transition tt^ at |g^| < 1 GeV^, it is shown in 

Fig. 2. 

With the photon wave function found, the decay form factors S* — > 77 and 
T — >■ 77 provide the opportunity to investigate the scalar and tensor meson wave 
functions. 



2.4.1 Form factor S l{q^)l 

Following the prescription given in previous sections for V l{<f)S and V 
^{q'^)P, we present the amplitude of the process S 7(g^)7 in terms of the light- 
cone variables. The transition form factor S reads: 

Fs^.rM\ 0) = ^^^^ j ^71^ / d'k^Ms)i^,is')Ss^,,is, s', q'). (37) 

Here, as before in Eq. (|2T|), the light-cone variables are introduced as follows: 
X = k2z/p, k± = k2±, and the gg invariant masses squared, s and s', are determined 
by Eq. (p^. The factor ■\/Nc, where Nc = 3 is the number of colours, is related to 
the normalization of the photon wave function performed in [^. 

The charge factor Zs-^-yy = 2^5^^^ is determined by the quark content of the 5- 
meson. We have two loop diagrams with quark lines drawn clockwise and anticlock- 
wise: the factor 2 in the determination of Zs^^-y stands for this doubling. For the /q- 
meson, one has two components with different charge factors (nn^-yy = (e^ -|- e^) / 
and Css^77 = e^, while for ao-meson Cao^77 — i^u ~ ^d)/v^- 

The spin structure factor Ss-,'yy{s,s',q'^) is fixed by the three-point quark loop 
trace for the amplitude of Fig. lb, with transverse polarized photons: 

- Sp[-f^^{k[ + m)7^^(A;i + m)(-A;2 + m)] = Ss^^^is, s', q^) gi^ . (38) 

Here 7^^ and 7^^ stand for photon vertices, 7q"'" = (7^"^^ 7^, while g-^^ is determined 
by formula (3) with the following substitution q ^ P — P' and g' P' . 
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One has 

Ss^-yyis, s', q^) = —2m 
Partial width, Ts^-yy, is determined as follows: 

msTs^yy = ^J d^ps;q,q')\A^,\^ = 7ia^\Fs^^y{0,0)\' . (40) 

Here ms is the scalar meson mass, the summation is carried out over outgoing 
photon polarizations, the photon identity factor, 1/2, is written explicitly, and the 
two-photon invariant phase space is equal to (i$2(P5; Q, q') = l/(167r). 

2.5 Two-photon tensor meson decay T ^ 77 

The decay amplitude T — > 7_L(g^)7 can be considered quite analogousy to the am- 
plitude of the two-photon scalar meson decay. The light-cone representation for the 
form factor F^^^{q^, g'^ = 0) with H = 0,2 reads: 

4"J,M 0) = / ^^r^ J d-k,Ms)i',(^')S?SrM, ,■') (41) 

^ ^ 

Here we use the same notation as in (^); the charge factors for the tensor and scalar 
mesons are equal to one another, Zt^^^ = Zs^^-y- 

The spin structure factors are fixed by the vertex for transition T qq; we 
denote this vertex as T^jy. One has: 

sZIp = Mli^CK + ^ha^h + m)T,Xk2 - m)] = (42) 

sZ^ip, <i)s^'\p ^ p '^ q') + s^;up^ q)s^'\p ^ p ^ q% 

where 7^^ and 7^-*- stand for photon vertices, 7^-*- = g^t^'ja', and g^^ is determined 
by Eq. (3) with the following substitution q P — P' and q' —>■ P' . The moment- 
operators S^^J^^f^{P, q) and Sj^J^^f^{P, q) work also in the intermediate state momentum 
space. Recall that q = P — P' and P^ = s, P'^ = s', (f = g^, while the momenta k[, 
ki and ^2 are mass-on- shell. 

The vertex T^jj{k) in its minimal form reads: 

2 , - 

T^u{k) = k^'j^ + k^-i^ - -g^^k , (43) 
where k = ki — k2 and g^yPv = 0. 



4m^ — s + s' 



IJ2 



Ass q 



2{s + s')q^ 



{s - s'f - g4 



(39) 
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Spin structure factors S^^^s, s', q^) and S'(^)(s, s', q^) are calculated by projecting 



on the moment-operators S^^\jP, q): 



^^^n^,^',g^) = f^f^^^^ (44) 

Explicit expressions for the spin structure factors S^^\s, s', q^) and S^'^'>{s, s', q^) are 
rather cumbersome, and we do not present them here. 

The qq (2"*"+) state can be constructed in two ways, namely, with the qq orbital 
momenta L = 1 and L = 3 (the ^P2qq and ^F2qq states). The vertex T^^, of Eq. (^Hj), 
corresponding to the dominant P-wave qq state, includes also certain admixture of 
the F-wave qq state. 

The vertex for the production of pure qq (L = 1) state reads: 

ri^='^ = k,T, + KT, - Igi^iTk), T, = ^i- ^J^^^ , (45) 

where the operator selects the spin-1 state for the qq (see |1^ ^ for detail). 
The {L = 3)-operator for the ^F2qq state is equal to: 

^i^='^ = kMTk) - ^ {gi^iTk) + r,k, + r^k,) . (46) 

For the qq wave function of tensor mesons, we use a parametrization similar to that 
for scalar mesons, see Eq. (^5|). The parameters Ct and br are determined by the 
tensor meson charge form factor at small g^, Eq. (pBf): the charge form factor is 
given by Eq. (^). The spin factor for the tensor meson, ST{s,s',q'^) is defined by 
the quark loop trace as follows: 

^Sp[r^,(A;i + mh^{k[ + m)T;,(^2 - m)] = 2P^^St{s, s', q^), (47) 

The operator T^,^ is written for the initial state transition T — > qq, Eq. (^3]), while 
T^j, describes the production of the outgoing tensor meson qq T that requires 
the following substitutions in (^): ki — > k[ and P ^ P'. The tensor meson charge 
form factor is averaged over polarizations that results in the factor 1/5 in (|47|). 

Partial width, Tt^jj, is determined as 

If 1 

mrTT^^^ = - J d<!){pT;q,q')- W'^,»fi\^ (48) 



4 2 
-vra 



\\F^%,{0,0)\' + \F^%,{OM 



Here rriT is the tensor meson mass, the summation is carried out over outgoing pho- 
ton polarizations, the photon identity factor, 1/2, is written explicitly; the averaging 
over tensor-meson polarizations results in the factor 1/5. 
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3 Results 



In this Section we present the results of calculations of partial widths for the radiative 
decays 0(1020) ^ 7/o(980), 0(1020) 77;, 77;', 771°, 7ao(980) and two-photon de- 
cays /o(980) ^ 77, ao(980) ^ 77, a2(1320) ^ 77, /2(1270) 77, /2(1525) ^ 77- 

3.1 ^(1020) 7/o(980): the decay amplitude and partial 
width 

Here we calculate the branching ratio for the decay 0(1020) — > 7/o(980) assuming 
the qq structure of /o(980). 

3.1.1 Wave functions of 0(1020) and /o(980) 

We write the wave functions of 0(1020) and /o(980) as follows: 

^0(5) = {nn sin (pv + ss cos (pv)^(j,is) , (49) 
*/o(980)(-s) = (nncosv? + sssinv^) V'/„(98o)(s) , 

assuming similar s-dependence for the nn and ss components. For ip^{s) and 
'0/o(98O)('5) the exponential parameterization is used, Eq. (^Sf ). The radius squared 
of nn component in the 0-meson is suggested to be approximately the same as that 
of the pion: R'^{nn) ~ 10.9 GeV~^, while the radius squared for the ss compo- 
nent, i?0(ss), appears to be slightly less, -R0(ss) ~ 9.3 GeV~^, that corresponds to 
bfp = 2.5 GeV~^. As to /o(980), we vary the radius of nn component in the interval 
6 GeV"^ < i?2^(98o)(nn) < 18 GeV'\ 

3.1.2 Partial width 0(1020) 7/o(980) 

The amplitude y4(^^^jQ(98o)(0) is the sum of two terms related to the nn and ss 
components: 

^<^-7/o (0) = COS ^ sin y?i/F{!^"]^^(ggo) (0) + sin ^ cos ^vF^Z\fo{980) (0) • (50) 

In our estimations we put cosipv ~ 0.99 and, correspondingly, | siny^y] ~ 0.14; for 
the /o(980) we vary the mixing angle in the interval 0° < \(f \ < 90°. 

The results of the calculation are shown in Figs. 3 and 4. In Fig. 3a the values 
^j!!i/o(98o)(0) and Aj^'5^^^(9gQ)(0) are plotted versus radius squared, R%(^9so)- 

In Fig. 4 one can see the value BR(0 — > 7/0 (980)) at various ip. Shaded areas 
correspond to the variation of (fv in the interval —8° < (fv < 8°; the lower and 
upper curves of the shaded area correspond to the destructive and constructive 
interferences of aJ!^^^^(98o)(0) and A^'J^^^(98o)(0), respectively. 
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The measurement of the /o(980) signal in the 77r'^'7r'' reaction (SND Collabora- 
tion) gives the branching ratio BR(0 7/o(980)) = (3.5 ± 0.3tli) x 10"^ in 



in the analysis of jir^Tr^ and 77r+7r channels (CMD Collaboration) it was found 
BR(0 7/o(980)) = (2.90 ± 0.21 ± 1.5) x 10"^ |n[; the averaged value is given 
in [H: BR(0 7/o(980)) = (3.4 ± 0.4) x 10"^ In our estimation of the per- 
missible interval for the mixing angle (p, we have used the averaged value given 



by [|T2[, with the inclusion of systematic errors of the order of those found in |jTT 
BR(0 ^ 7/o(980)) = (3.4 ± 0.41^;^) x 10"^ 

The calculated values of BR(0 7/0 (980)) agree with experimental data for 
\(f \ > 25°; larger values of mixing angle, \ > 55°, correspond to a more compact 
structure of /o(980), namely, -^^^(gso) ^ 10 GeV~^, while small mixing angles \(p\ ~ 
25° are related to a loosely bound structure of the /o(980), -^^^(gso) > 12 GeV~^. 



The evaluation of the radius of /o(980) was performed in |19] on the basis of 



GAMS data where the t-dependence was measured in the process 7i~p — > 
/o(980)n (t is the momentum squared transferred to /o(980)): these data favour a 
comparatively compact structure of the qq component in /o(980), namely, -^^^(ggo) = 
6 ± 6 GeV-2. 



3.2 Radiative decays 0(1020) -fi], 77/', 77r^ 7ao(980) 

The decays 0(1020) — 777,777' do not provide us with a direct information on the 
quark content of /o(980) and 0(1020); still, calculations and comparison with data 
are necessary to check the reliability of the method. The decays 0(1020) — > 771° and 
0(1020) — >■ 7ao(980) allow us to evaluate the admixture of the nn component in the 
meson; as is seen in the previous section, this admixture affects significantly the 

value r0(io2O)-^7/o(98O)- 

For the transitions — > 777 and — > 777' we take into account the dominant ss 
component only: — sin^ss in 77-meson and cos6'ss in r/'-meson, with sinO = 0.6. 

For the pion wave function we have chosen 6^ = 2.0 GeV^^ that corresponds to 
R^. = 10.1 GeV~^, the same radius is fixed for the nn component in rj and t]'. As to 
the strange component in r] and r]', we put its slope to be the same: br,{ss) = br^'(ss) = 2 
GeV~^, that leads to a smaller radius i?^(ss) = 8.3 GeV~^. 

The calculation results for branching ratios compared to those given by PDG- 
compilation |jl2| follows: 

BR(0 r/7) = 1.46 x 10"^ , BRpdg(0 Vl) = (1-30 ± 0.03) x 10^^ ^ 
BR(0 77'7) = 0.97 X 10-^ , BRpdg(0 v'l) = (0.67l°f ) x 10^^ (51) 
As is clearly seen, the calculated branching ratios agree reasonably with those given 

inig. 

For the process 77r° the compilation [|T^ gives BR(0 77?*^) = (1.26 ± 
±0.10) X 10"'^, and this value requires | sinyjyl ^ 0.07 (or \ipv\ — 4°), for just with 



15 



this admixture of the nn component in 0(1020) we reach the agreement with data. 
However, in the estimation of the allowed regions for mixing angle ip, Fig. 4, we use 

|(^y|=4°±4° (52) 

considering the accuracy inherent to the quark model to be comparable with the 
obtained small value of |v?y|- 

The process 0(1020) 7ao(980) depends also on the mixing angle |v?y|: the 
decay amplitude is proportional to siny^y, namely, A^_,^ao = sin (/^yA^!^^^^. For the 
region -Rao(980) ~ ^ GeV~^ - 12 GeV~^, our calculation gives the following branching 
ratio: 

BR (0(1020) ^ 7ao(980)) = sin^ ipy ■ (14 ± 3) x 10^^ (53) 

with lower values for R^aQ{9m) ~ ^ GeV"^ and larger ones for -^^^(gso) ~ 12 GeV~^. 
At sin^ = 0.01 ± 0.01, we have BR(0(1O2O) ^ 7ao(980)) = (0.14 ± 0.14) x 10"^ 



In |11], the rjir^ spectrum was measured in the radiative decay 0(1020) — > '-yrjTt^: 
it was found that BR(0(1O2O) ^ 77/71°; M^^ > 900 MeV) = (0.46 ± 0.13) x 10"^ 



This value does not contradict the equation ( ^31) with \>fv\ = 8°; moreover, if the 
ratio of the background/resonance in the region M^^r ~ 900 MeV is not small, that 



is rather possible, the value found in |11] agrees with smaller values of Iv^vl- 



3.3 Radiative decays /o(980) 77 and ao(980) 77 

The amplitude (980)^77 (0, 0) is determined by contributions of two flavour com- 
ponents: 

^/o(980)^77(0, 0) = COS 0)^77(0' 0) + ^Ffom)^yyi^' 0) ' (^4) 

The amplitudes ^^"(930)^77(0, 0) and Aj^(ggQ-)^^^(0, 0) depend on the radius squared 
of /o(980): these amplitudes plotted versus -^^^(ggo) are shown in Fig. 36. 

Figure 5 demonstrates the comparison of calculated partial width r/Q(98o)^77, at 
different -R^q(98o) ^.nd with the magnitude found in [l^: rjo(98o)^77 = 0.28to;?3 . 
It is possible to describe the data using positive mixing angles 77° < (p < 93° as 
well as negative ones: (—54°) < (p < (—38°). 

The amplitude for the decay ao(980) — * 77 is determined by the similar form 
factor that is similar to that of nn component in the /o(980), with the only difference 
C/o — >■ Cao- The amplitude ^4^^ (980)^77 /2Cao{980)->77 is shown in Fig. 36 as a function 
of -Rap(98o)- Experimental study of r(ao(980) 77) was carried out in Refs. |21, 
13, the averaged value is: T{r]7r)T{^^)/Ttotai = 0.24;°;°? keV 0. Using Tt^tai ^ 
T{f]7i) + T{KK), we have r(ao(980) ^ 77) = 0.30t{j:}J keV. The calculated value of 
r(ao(980) 77) agrees with data at Rag(980) belonging to the interval 10 GeV~^ < 
^ao(98o) < 26 GeV-2: the values of Rl^^^gso) the order of ~ 10 - 17 GeV'^ look 
quite reasonable for a meson of the l^Poqq multiplet. 
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3.4 Partial widths for the two-photon decays of tensor mesons 



Here we present our results for the decays 02(1320) — > 77, /2(1270) — > 77 and 
/2(1525)-.77- 



3.4.1 Decay 02(1320) 77 

a2(1320)->77^ 



( H) 

The form factor fi320W'Y'Y('-'' 0) is equal to that for the nn component, up to the 



charge factor: 



-^if(1320)^77(0' 0)/^a2(1320) - -F^fl^^ (0, 0) , (55) 



This universal form factor is determined by the spin structure of the vertex 02(1320) — 
qq, which is regulated by the admixture of the ^F2gg-state in O2(1320). Figure 6 



demonstrates the calculated form factor Fn^l^^^{0,0)/Znn as a functions of for 
different vertices given by Eqs. (|43|), (^5]) and (|^). 

The comparison with data is performed for the form factors calculated with 
the use of minimal vertex given by Eq. (P3|). In Fig. 7 we plot calculated values of 
ra2{i320)^77 being a function of -^^2(1320) versus the data: ra2(i32o)^77 = 0.98±0.05± 
0.09 keV and 1^3(1320)^77 = 0.96 ± 0.03 ± 0.13 keV The calculated value 

of 

(1320)^77 reproduces data with Rl^(^i^20) ~ ^ GeV ^ only. Still, one should 
not be convinced that larger values of the radius -^^2(1320) excluded by the data. 
Experimental extraction of the signal O2(1320) 77 faces the problem of a correct 
account for coherent background. This problem has been investigated in [^: it was 



shown that the measured value of (1320)^77 can fall down by a factor ~ 1.5 due 
to the interference "signal-background". Therefore, being careful, we estimate the 
region for ra2(i32o)^77 allowed by the data as 1.12 keV< ra2(i32o)^77 < 0.60 keV. 

Comparing the calculations with data, one should take into account uncertainties 
inherent to the model. In our calculation, ra2(i32o)-^77 strongly depends on the 
constituent quark mass. In Fig. 6, the form factors and partial widths are depicted 
for rriu^d = 350 MeV and = 500 MeV. However, decreasing of constituent quark 
mass by 10% results in the increase of the form factor F^^^(0,0) approximately by 
10%, that means the 20% growth of the calculated value of ra2(i32o)— >77 at fixed 
-^02(1320)- '^^^ uncertainty in the definition of the constituent quark mass looks 
quite reasonable, therefore, the 20% error in the model prediction for (1320)^77 is 
to be regarded as normal. 

Summing up, the calculation of Fa2(i32o)^77 with the minimal vertex O2(1320) 
qq given by (^) provides reasonable agreement with data at 7 GeV~^ < -Ra2(i32o) ~ 
13 GeV-2. 



The vertex corresponding to the production of qq pair in the F-wave, Eq. (^6]), 
gives the partial width value ~ 0.1 keV that is by an order of value less than for the 
P-wave qq component. 
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The vertex for pure P-wave qq state, see (^5]), gives us a 10% smaller value of 
ra2{i320)^77 ^s Compared to what is provided by minimal vertex (^31). 



3.4.2 The decays /2(1270) 77 and /2(1525) ^ 77 

We define the wave functions of /2(1270) and /2(1525) as follows: 

^/2{i270)(s) = (cosv?T nn + siny^T ss) ^t(s), 

^/2(i525)(s) = {- sin ipT nn + COS ipT ss)ijT{s). (56) 

Then, the form factors for the two-photon decays of /2-mesons read: 

^hL70)->77(0' 0) = COS y^TFinl^-y{0, 0) + sin ^TF,f2^^{0, 0), 
^/S525H77(0,0) = -smipTFi^l,,{0,0) + cos^TF^^2,,{0,0). (57) 



Following [0, IJ, |25[ , we put the following values for partial widths: T j2 (1270)^77 = 
(2.60 ± 0.25t[]:^^) keV and 1/3(1525)^^,^, = (0.097 ± 0.015t[}:^g) keV. The magnitude 
of the extracted signal depends on the type of model used for the description of 
the background. With coherent background, the magnitude of the signal decreases, 
and the second error in F/j (1270)^77 and F/2(i525)^^^ is related to the background 
uncertainties. 

A satisfactory description of data has been reached with R"^ < 10 GeV~^ and 
two mixing angles ipT'- — 0° and (fx — 25°, see Fig. 8. For example, at 
-Rjj, = 9 GeV~^ and ipT = 0°, the calculations result in Fj2(i27o)^77 = 2.240 keV and 
r/2{i525)-*77 = 0.090 keV. Nearly the same values are reproduced at Rj^ = 9 GeV~^ 
and = 25°, namely, Fj2(i270)^77 = 2.237 keV and Fj2(i525)^77 = 0.093 keV. 

As for the reaction 02(1320) — > 77, the model uncertainties, ~ 20%, are inherent 
in calculations of Fj2(i27o)-^77 and Fj^ (1525)^77- 



4 Conclusion 

Figure 9 demonstrates the (y?, i?j^(ggQ-))-plot where the allowed areas for the reac- 
tions 0(1020) — *• 7/0 (980) and /o(980) — > 77 are shown. We see that radiative 
decays (^(1020) — > 7/o(980) and /o(980) 77 are well described in the framework 
of the hypothesis of the dominant qq structure of /o(980). The solution with neg- 
ative (p seems more preferable. For this solution the mixing angle cp for nn and ss 
components (rin cos y^ + ss sin </)) is equal to (f = —48° ±6°, that is, the qq component 
is rather close to the flavour octet (v^octet = —54.7°). However, the radiative-decay 
data do not exclude the variant when the /o(980) is almost pure ss state with 
(f = 85° ±5°. 

The dominance of the quark- ant iquark component does not exclude the existence 
of other components in /o(980) on the level 10% — 20%. The location of resonance 
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pole near the KK threshold definitely points to certain admixture of the long-range 
KK component in /o(980). To investigate this admixture the precise measurements 
of the KK spectra in the interval 1000 — 1150 MeV are necessary: only these spectra 
could shed the light on the role of the long-range KK component in /o(980). 

The existence of the long-range KK component or that of gluonium in the 
/o(980) results in a decrease of the ss fraction in the qq component: for example, if 
the long-range KK (or gluonium) admixture is of the order of 15%, the data require 
either = -45° ± 6° or v? = 83° ± 4°. 

There is no problem with description of the decay ao(980) 77 within the 
hypothesis about the qq origin of the ao(980): the data are in a good agreement 
with the results of calculations at -^^^(gso) ~ 10 - 17 GeV~^. 

We have calculated the two-photon decays of tensor mesons, members of the qq 
multiplet l^P2qq. The calculated partial widths of radiative decays, 02(1320) 77, 
/2(1270) 'j'j and /2(1525) 77, are in a reasonable agreement with data. The 
radial wave functions of 02(1320), /2(1270) and /2(1525) are close to those of ao(980) 
and /o(980) found in the study of radiative decays ao(980) — > 77, /o(980) 77 and 
0(1020) — > 7/o(980). The possibility to describe simultaneously scalar and tensor 
mesons using approximately equal radial wave functions may be considered as a 
strong argument that all these mesons, tensor ones, a2(1320), /2(1270) and /2(1525) 
and scalar ones, ao(980) and /o(980), are members of the same P-wave qq multiplet. 

We are grateful to L.G. Dakhno, V.N. Markov, M.A. Matveev and A.V. Sarantsev 
for useful discussions. The paper is suppored by RFBR grant 01-02-17861. 
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Figure 1: a) Diagrammatic representation of the transition 0(1020) — > 7/o(980). b) 
Three-point quark diagram: dashed hues I and II mark two cuttings in the double 
spectral representation. 




Figure 2: Photon wave function for non-strange quarks, ip^^nn{k'^) = 9'y{k'^)/{k'^ + 
m^), where = s/4 — m^; the wave function for the ss component is equal to 
ip^^ssik"^) = g'y{k^)/{k^ + "^s); the constituent quark masses are m=350 MeV and 
m^^SOO MeV. 
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Figure 3: Amplitudes for strange and non-strange components, ss and nn, 
as functions of the /o(980)-meson radius squared: a) F^}^f^{0)/Z^}^j!^ and 
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Rfo(98o) (nn), GeV" 



Figure 4: Branching ratio BR{(f){1020) 7/o(980)) as a function of radius squared 
of the nn component in /o(980). The band with vertical shading stands for the 
experimental magnitude: BR(0 ^ 7/o(980)) = (3.4 ± OAtH) x IQ-^ Five other 
bands, with skew shading, correspond to |(/?| = 24°, 37°, 53°, 64°, 80° at -8° < cpv < 
8°. 
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(nn), GeV" 




R'o(98o) (nn), GeV" 



Figure 5: Partial width rjj,(98o)-+77; experimental data are from [15] (shaded area), 
a) Curves are calculated for positive mixing angles </? = 77°, 80°, 85°, 90°, 93° and b) 
negative angles ip = -38°, -43°, -47°, -54°. 
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Figure 6: Transition form factors T — > 77 (see (17) or (24)) for the non-strange (nn) 
and strange (ss) quarks versus mean tensor meson radius squared, R^. a) -^^^^(O, 0) 
and fIi^\o,0) for l^P2qq state with minimal vertex, Eq. (41); b) the same as Fig. 
6a but with the vertex determined by (45); c) Fgq\o,0) and Fgq\o,0) for l^F2qq 
state with vertex given by (46). 




Figure 7: Partial widths for a2 — > 77 versus mean tensor-meson radius squared, 
B^. Thick sohd hne: r(a2(1320) '-ff) for the vertex given by (43); Dashed hne: 
1(02(1320) 77) for the vertex given by (45) Dotted hne: T{a2{^ 2000) 77) for 
the vertex given by (46). 
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Figure 8: The {R^, </')-plot, where ip is mixing angle for flavour components 

/2(1270) = nncosif + ss sin (/? and /2(1525) = —nnsimp + sscos(/9, with hatched 
areas which show the regions allowed by data for decays /2(1270) 77 and 
/2(1525) ^ 77. 
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Figure 9: The ((/?, R'j^^(^ggQ^)-p\ot: the shaded areas are the allowed ones for the reac- 
tions 0(1020) ^ 7/o(980) and /o(980) ^ 77. 
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